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Abstract 
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Climate change can exacerbate future regional air pollution events by making condi- 
tions more favorable to form high levels of ozone. In this study, we use spectral nudging 
with WRF to downscale NASA earth system GISS modelE2 results during the years 
5 2006 to 2010 and 2048 to 2052 over the continental United States in order to com- 
pare the resulting meteorological fields from the air quality perspective during the four 
seasons of five-year historic and future climatological periods. GISS results are used 
as initial and boundary conditions by the WRF RCM to produce hourly meteorological 
fields. The downscaling technique and choice of physics parameterizations used are 
io evaluated by comparing them with in situ observations. This study investigates changes 
of similar regional climate conditions down to a 12 km by 12 km resolution, as well as 
the effect of evolving climate conditions on the air quality at major US cities. The high 
resolution simulations produce somewhat different results than the coarse resolution 
simulations in some regions. Also, through the analysis of the meteorological variables 
15 that most strongly influence air quality, we find consistent changes in regional climate 
that would enhance ozone levels in four regions of the US during fall (Western US, 
Texas, Northeastern, and Southeastern US), one region during summer (Texas), and 
one region where changes potentially would lead to better air quality during spring 
(northeast). We also find that daily peak temperatures tend to increase in most major 
20 cities in the US which would increase the risk of health problems associated with heat 
stress. Future work will address a more comprehensive assessment of emissions and 
chemistry involved in the formation and removal of air pollutants. 
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1 Introduction 
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Changes in climate, emissions, population, technologies, and land-use can impact air 
25 quality in a variety of ways. Studies suggest that climate change can exacerbate future 
regional air pollution events (e.g. (Liao et al., 2006; Mickley et al., 2004; Stevenson et 
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al., 2006; Weaver et al., 2009) by making conditions more favorable to form high lev- 
els of ozone, e.g. by increasing temperature and biogenic emissions and decreasing 
ventilation. Increased temperatures affect air quality by affecting reaction rates and gas 
solubility in water droplets (EPA, 1989). Pollutant dispersion and removal are affected 
5 by large-scale circulation patterns and precipitation, while cloud cover frequency and 
duration impacts photolytic activity, which in turn affects reaction rate coefficients and 
conversion of gases to aerosols. Stagnation event frequency and duration affects the 
mixing of polluted air with air above the boundary layer. To simulate regional air quality 
dynamics, results from global models are downscaled using dynamical downscaling. 
io The starting point of dynamical downscaling is typically a set of coarse- resolution large- 
scale meteorological fields (either from a General Circulation Model, GCM, or from 
global reanalysis data) which are used to provide the initial, and lateral and surface 
boundary conditions to a regional climate model (RCM). Typically the RCM simulation 
does not feedback into the GCM, but adds regional detail in response to finer scale 
15 forcing (e.g. topography, land use/land cover) as it interacts with the larger-scale atmo- 
spheric circulation (Giorgi, 2006). In this study, we address the benefits of downscaling 
using an RCM when analyzing the implications of climate change on air quality and 
health, especially in urban areas. 

Recently, climate modeling efforts have shifted their focus from analyzing mean val- 
20 ues of climate variables (e.g. temperature and precipitation) to extreme values, variabil- 
ity, and shifts in the frequency of climate patterns that are more relevant for air quality. 
Jacob and Winner (2009) compiled results from a number of studies on climate change 
and air pollution, and summarized that increases in regional stagnation consistently in- 
creases ozone and particulate matter (PM) concentrations. They show that positive 
25 temperature perturbations consistently increase ozone while they can sometimes de- 
crease PM concentrations and that positive perturbations in mixing depth, wind speed, 
cloud cover, and precipitation all decrease ozone and PM concentrations to varying 
degrees. Mickley et al. (2004a) applied the Goddard Institute for Space Studies (GISS) 
GCM 2’ (Rind and Lerner, 1996; Rind et al., 1999) with implemented carbon monoxide 
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(CO) and black carbon (BC) tracers to simulate the impact of climate change on air 
quality in the US. They found that increased severity of future pollutant episodes in 
the northeast and midwest during the summer was due to a decrease in frequency of 
surface cyclone tracking in southern Canada. 

Dawson et al. (2008) developed a global-regional climate-air pollution modeling sys- 
tem (GRE-CAPS) by coupling GISS II GCM, MM5 regional meteorological model, and 
the PMCAMx regional CTM. This system was intended to enable studies of the effects 
of changes in climate, intercontinental transport, and emissions on regional and urban 
air quality. Leung and Gustafson (2005) applied a similar approach to assess the poten- 
tial effects of climate change in the United States. They developed meteorological fields 
by downscaling the NASA GISS GCM simulations using an MM5-based RCM (Grell, 
1994). Their analyses were based on changes in surface air temperature and down- 
ward solar radiation, precipitation frequency, stagnation events, and ventilation. They 
defined a time to be stagnant when for four consecutive days the following criteria are 
met: (a) the 10m wind speed is less than 4 ms -1 (b) the 500 mb wind speed is less 
than 13 ms -1 at 07:00 LST, and (c) the total rainfall is less than 0.001 cm for the 4 day 
period (Korshover and Angell, 1982). They also compared the daily average number of 
unvented hours, which are hours when the product of the mean wind speed within the 
boundary layer and the boundary layer height is less than 6000m 2 s _1 (Pielke et al., 
1991). In Tagaris et al. (2007), meteorological inputs to the Community Multi-scale Air 
Quality (CMAQ) Chemistry/Transport Model (CTM) were developed using the meteo- 
rological fields of Leung and Gustafson (2005) to investigate the potential impacts of 
global climate change and emissions on regional air quality using CMAQ (Byun and 
Schere, 2006). Similarly, Nolte et al. (2008) investigated the impact of climate change 
on future air quality in the United States by dynamically downscaling outputs from the 
GISS GCM with the MM5 RCM and predicted an increase in 0 3 over Texas and large 
portions of the southeast using CMAQ model. 

These studies have illustrated the value of using the regional downscaling approach 
in order to better understand the impact of climate change on regional air quality. With 
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the strong dependence on localized flow patterns, air quality models benefit from the 
higher-resolution wind, temperature, precipitation, and boundary layer structures pro- 
duced by a RCM (Leung and Gustafson, 2005). Weaver, et al. (2009) stresses that the 
science of coupling global climate and regional air quality models is still at a young 
5 state and that there are particular questions as to which climate metrics and statistics 
are most relevant to air quality and how sensitive simulation results are to downscaling 
methodologies. 

In our previous recent work (Liu et al., 2012) we examined the performance of two 
nudging techniques, grid and spectral nudging, by downscaling NCEP/NCAR data us- 
io ing the Weather Research and Forecasting (WRF) Model and identified benefits of 
spectral nudging at producing small scale features while preserving the large scale 
forcings. Following these findings, in this study, we use spectral nudging to downscale 
the NASA earth system GISS modelE2 results during the years 2006 to 2010 and 
2048 to 2052 over the continental United States (CONUS) in order to compare the re- 
15 suiting meteorological fields from the air quality perspective during the four seasons of 
five year historic and future climatological periods. GISS results are used as initial and 
boundary conditions by the WRF RCM to produce hourly meteorological fields. The 
downscaling technique and choice of physics parameterizations used are evaluated by 
comparing them with in situ observations. This study investigates changes of similar 
20 regional climate conditions down to a 12 km by 12 km resolution, as well as the effect 
of evolving climate conditions on the air quality at major US cities. 

2 Approach 

(/) 

In this study a regional climate model is used to downscale a global climate model to 
develop meteorological fields for the present and future. Each component of the mod- 
25 eling system is described below along with the downscaling and evaluation methods 
used. 
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2.1 Model descriptions 

O 

Global model 

o 

Lateral boundary and initial conditions for the regional forecast modeling are taken 
from the GISS ModelE2. The model has a horizontal resolution of 2° x 2.5° latitude by 
5 longitude. The vertical discretization has 40 layers and follows a sigma coordinate up 
to 150 hPa, with constant pressure layers between 150 and 0.1 hPa. The surface is 
split into four types: open water (including lakes, rivers and oceans), ice-covered water 
(including lake ice and sea ice), ground (including bare soil and vegetated regions) and 
glaciers. 

io Simulations are carried out for the calendar years 2006-201 0 and 2048-2052, driven 
by future atmospheric conditions over the 21st century and follow the scenario devel- 
opment process for IPCC AR5. The specific scenario used for this study is the “Rep- 
resentative Concentration Pathway” (RCP) 4.5 (Lamarque et al., 2011; Moss et al., 
2010), that is a scenario of decadal global emissions of greenhouse gases, short- 
15 lived species, and land-use-land-cover which produces an anthropogenic radiative 
forcing at 4.5Wm' 2 (approximately 650 ppm C0 2 -equivalent) in the year 2100. The 
detailed characteristics of this scenario are enumerated in Moss et al. (2010). The at- 
mosphere/terrestrial biosphere-only version of the GISS modelE2 was driven by sea- 
ice and sea-surface temperature conditions calculated by the coupled earth system 
20 model version that is submitted to the CMIP5 archive. The model spinup time was 3 
years, starting either from 2003 or 2045, which is sufficient for the dynamic equilibration 
of the model’s climate and chemically active tracers. SST and sea-ice boundary con- 
ditions vary both seasonally and interannually, GHG concentrations change annually, 
and emissions change annually by linearly interpolating the decadal CMIP5 emission 
25 datasets. 6-hourly instantaneous outputs of physical and chemical parameters were 
produced for regional downscaling by WRF (Sect. 2.2). 
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Regional model 

The regional climate model used is the non-hydrostatic Weather Research and Fore- 
casting (WRF) Model (Skamarock and Klemp, 2008) version 3.4. The simulation do- 
main covers the CONUS and portions of southern Canada and northern Mexico and 
is centered at 40° N and 97° W with dimensions of 164 x 138 horizontal grids cells with 
a grid-spacing of 36 km. It contains 35 vertical levels, with the top pressure of 50 hPa. 
The configuration of physics schemes is as follows: the long-wave Rapid Radiative 
Transfer Model (RRTM) (Mlawer et al., 1997) and Dudhia scheme (Dudhia, 1989) are 
used for longwave and shortwave radiation respectively; the Yonsei University (YSU) 
(Hong et al., 2006) scheme is used for the planetary boundary layer; the Noah scheme 
(Ek et al., 2003) is used for land surface model (LSM); a revised version of Kain-Fritsch 
scheme (Kain, 1993) is used to represent the effects of both deep and shallow cumulus 
clouds; Lin et al. (1983) is chosen for cloud microphysics. 

2.2 Dynamical downscale of global results 

The GISS ModelE2 fields include temperature, relative humidity, horizontal wind veloc- 
ities, soil temperature and moisture at different soil depths, sea surface temperature, 
surface pressure, ice fraction and snow water equivalent. The WRF Preprocessing Sys- 
tem (WPS), which reads in this global data and interpolates it to the WRF grid points, 
does not process GISS data directly. Therefore, an interface program was developed 
to link the GISS output with WPS. 3-D variables, such as wind and temperature, are 
interpolated to 21 fixed pressure levels; the lowest level of these 3-D variables and 
surface level properties were vertically interpolated to produce 2 m temperature, 2 m 
humidity, and 10 m wind fields. The soil-related variables were also interpolated to the 
depths defined from the LSM. 

Global model results are used as initial and boundary conditions for the regional cli- 
mate simulations. Spectral nudging with a wave number of 3 in both zonal and merid- 
ional directions is used, i.e. all waves with wave numbers greater than 3 are filtered (Liu 
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et al., 2012); no nudging is conducted for shorter waves to provide similarity with the 
large scale GCM simulation but allow small scale features to freely develop (Liu et al., 
2012). Spectral nudging is applied to temperature, horizontal winds, and geopotential 
heighst. Only horizontal winds are nudged at all vertical levels, while no nudging is 
5 conducted for other variables within the planetary boundary layer (PBL). The nudging 
coefficient for all variables was set to 3 x 1 0 -4 s -1 (Stauffer and Seaman, 1 990). During 
the simulation, nudging is conducted every 6h, consistent with the frequency of the 
GISS data. 

D 

c ri 
o 
c: 

3 Model application and evaluation 

13 

io WRF is applied here using a nested grid approach. The modeling domain uses a Lam- 
bert Conformal Projection centered at 40° N, 97° W with true latitudes of 33° N and 
45° N. The outer domain uses a 36 km horizontal grid-spacing that covers the entire 
continental US as well as portions of Canada and Mexico (5940 x 5004 km). Two inner- 
most domains cover 984 x 1020km , and, 948 x 948km regions with 12 km horizontal 
is grid-spacing and focusing on the northeast and southeast US respectively (Fig. 1 ). The 
periods modeled are 2006 through 2010 (historic) and 2048 to 2052 (future). The sim- 
ulated coarse-grid hourly meteorology is used as initial and boundary conditions for the 
finer grids. 

Observations are used to evaluate the ability of GISS-WRF to reproduce the long- 
20 term yearly climatic means, and the meteorological fields that strongly impact air 
quality. The model performance is evaluated by using statistical measures. This is 
a common analysis that is proposed by Emery and Tai Emery (2001) and has been 
adopted by the meteorological modeling community. Statistics such as mean bias (MB), 
mean absolute gross error (IMAGE_Team) and root mean square error (RMSE) are 
25 calculated: 
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5 where P; is the predicted value of the tested parameter (i.e. temperature), O, is the 
corresponding observed value, and N is the total number of the predictions used for 
the comparison. MAGE gives an estimation of the overall discrepancy between pre- 
dictions and observations, while MB is sensitive to systematic errors. The root mean 
square error (RMSE) incorporates both the variance of the prediction and its bias. 
io Additional details for the above evaluation metrics can be found in Yu et al. (2006) 
The observations used for the statistical analysis are TDL (Techniques Development 
Laboratory) data from the Research Data Archive (RDA) http://dss.ucar.edu in dataset 
number ds472.0), which is maintained by the Computational and Information Systems 
Laboratory (CISL) at the National Center for Atmospheric Research (NCAR). These 
is are hourly surface observations for wind speed, wind direction, and temperature dur- 
ing the four seasons over a five year period (2006-2010). In the statistical analysis, 
the continental US domain has been divided in 4 sub regions, the western (W), the 
midwest (MW), the south (S) and the northeast (NE) USA (Fig. 1 ). 


4 Results 

20 Table 1 summarizes the comparison of the GISS-WRF modeling system predictions 
for wind speed and direction against TDL hourly surface observations during the four 
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seasons of a five year period (2006-2010) over four regions of the continental USA: 
Western USA, the Midwest, the South, and the Northeast. Overall the model predic- 
tions agree well with observations with the MB over the total domain ranging between 
-0.1 ms -1 (during summer) 0.2 ms -1 (during spring). The model performance is better 
during summer with RMSE as a low as 2.2 ms -1 over the south and worst during winter 
with RMSE over the West up to 3.9 ms -1 . Wind speed is better predicted over south 
(with MAGE ranging from 1.7 ms -1 to 2.2 ms -1 ) while wind direction is better predicted 
over northeast (with MAGE ranging from 72 deg. to 78 deg.). Table 2 summarizes the 
comparison of the GISS-WRF modeling system predictions for temperature against 
observational data during the four seasons over the continental USA. Compared to 
observations, the model tends to under predict temperature during winter (MB up to 
-7.5 K), spring (MB up to -2.7 K), and summer (MB up to -1.9 K over West) but over- 
predict temperature during fall (MB up to 2.9 K). These biases are mostly over the 
Western US and correspond to the biases in the GISS fields. Model performance is 
better over the Southern US, especially during summer (RMSE = 3.5 K). 

D 

4.1 Temperature 

C 

O) 

The 5 yr mean of the modeled 2 m air temperature across the simulation domain for the 
future is 1 K warmer than that of the historical simulation (284 and 285 K respectively) 

(Fig. 2). Consistent with other studies (Leung and Gustafson, 2005; Liao et al., 2007; 

Nolte et al., 2008; Tagaris et al., 2007; Woo et al., 2008) most of the warming, between 
3 and 4K, occurs over the western states (California, Nevada, Arizona, Texas, and 
Utah) and over western Canada (Fig. 2a) and the results of a t test suggest that the 
warming in this region is statistically significant (p value < 0.05 in this region). Signifi- 
cant warming mainly occurs over these regions during the winter and spring months, 
where average surface temperature change reaches 4 degrees in western states, es- 
pecially in and around Nevada (the p value is less than 0.05 for these temperature 
increases). Since temperatures are low during the winter and spring, warming during 
these seasons may not lead to increased concentrations of secondary pollutants such 

2526 
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as 0 3 and secondary PM, but warming could lead to decreased emissions of PM from 
heating processes such as wood burning (e.g. from wood stoves). The GCM simu- 
lations predict a similar warming pattern during the winter and spring, but only up to 
3 degrees K (Fig. SI). During the summer months, Texas and northeastern Canada 
5 experience a warming of 2 and 4 degrees respectively, although the GCM predicts up 
to one degree more warming over western Texas and the p value associated with the 
downscaled temperature changes over western Texas is between 0.05 and 0.10. An 
average warming of 3 degrees occurs over the Midwest (p value < 0.05) and a warm- 
ing of around 2 degrees also occurs over most of Texas and Eastern US (the p value is 
10 between 0.05 and 0.15 which is not significant) during the fall. The eastern US states, 
on the other hand, are cooler during the winter and spring months with the southeast- 
ern states and Texas cooling up to just less than 2 degrees, however the cooling here 
is not statistically significant (p value greather than 0.05). 

The smaller, more highly resolved, 12 km domain over the northeast simulates sim- 
15 ilar magnitudes of temperature change to the 36 km domain. The root mean square 
difference of the future temperature change between the 36 km and 12 km domains is 
very small (less than 0.004 K); indicating the similarity between the two simulations. 

The standard deviations of the simulated temperature changes in the northeast for the 
12 km (standard deviation of 0.25 K 2 ) and 36 km (0.24 K 2 ) domains show that the fine 
20 resolution simulation introduces slightly more variability than the coarse resolution do- 
main, especially during the winter (0.49 K 2 for the 36 km and 0.52 K 2 for 12 km) and 
spring (0.27 K 2 for 36 km and 0.29 K 2 for 12 km). The northeast sees cooling of less 
than 1 degree during the spring and warming of up to 2 degrees during the summer 
(Fig. 2b). During fall, large warming between 2-3 degrees is simulated over much of 
25 New York State. Over the southeastern 12 km domain, similar warming occurs as the 
36 km domain, ranging between 1 and 3 degrees during the summer and fall with the 
greatest warming occurring during the fall over North Carolina and Tennessee (Fig. 2c). 
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4.2 Insolation and precipitation 

O 

C 

A change in downward solar radiation at the surface, or insolation, is an indicator of 
changes in cloudiness. For this reason, spatial distributions of the change in insolation 
at the surface are similar in structure to average daily precipitation, but not identical. 

5 Spatial distributions of surface temperature and insolation changes have similar struc- 
tures in some cases. Weaver et al. (2009) explain that these meteorological conditions 
can have either competing or reinforcing effects on air quality. When temperature and 
insolation change in the same direction, 0 3 concentrations tend to change in the same 
corresponding direction, whereas temperature and insolation varying in opposite direc- 
io tions correspond with mixed changes in 0 3 . 

Decreases in daily mean precipitation are found over the Pacific coast where some 
regions receive 2 mm less rain per day (or 30% less rain), on average, and some de- 
creases were simulated over the southeastern region (Fig. 3a). Reduced rain along the 
Pacific coast occurs mostly during the winter, as a major portion of Western US sees 
is greater than a 2 mm per day decrease. Correspondingly, insolation over the Pacific 
coast increases during the winter by up to 15Wm -2 (Fig. 4a). The southeast experi- 
ences a similar magnitude of drying, but mainly during the fall. Both the 36 and the 
12 km simulations over the southeast predicted greater than 2 mm less rain per day 
during the fall (Figs. 3c and 4c), which is also consistent with insolation changes in the 
20 region (increase of up to lOWrrf 2 ). Interestingly, the high-resolution simulation pre- 
dicts that the southeast receives up to 2 mm per day more rain during the summer, 
which is not apparent in the 36 km domain. The 36 and 12 km resolution simulations 
over the northeast on the other hand, predict more rain over most of Vermont, New 
Hampshire and Maine during the summer, while most of Connecticut and New York 
25 receive less rain (Fig. 3b). The precipitation trend in the northeast reverses during the 
fall when Connecticut and New York receive more rain and the states farther north are 
dryer. There is a decrease in insolation of 5-15 Wm' 2 during the spring and fall and an 
increase during the fall of 10-20 Wm -2 over the northeast (Fig. 4b). Correspondingly, 
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the temperature only decreases slightly during spring, while it increases during sum- 
mer. However, contrary to the insolation trend, the temperature in the northeast sees 
large increases during fall. 

=3 

4.3 Stagnation events 

CD 

— S 

5 Stagnation events occur when wind speeds are low and little precipitation occurs over 
an extended period. Since transport and deposition of pollutants is decreased during 
a stagnation period these events promote poor air quality. During the winter and spring 
months, the spatial distribution of the number of stagnation days per season does 
not change significantly over the US (Fig. 5a). Over southern Texas, the number of 
io stagnation days during fall increases in some small areas by 1 0 to 1 5 days per season, 
which correlates with the increase in temperature in the region. Large portions of this 
region already see over 30 days of stagnation per season. Stagnation decreases over 
Texas during the summer. Over most of the southeast stagnation days also decrease by 
up to 10 days per season corresponding to the increase in precipitation (Fig. 5c) which 
is is large compared to the average number of stagnation days during the summer of the 
historic simulation (between 1 5 and 30). While the 36 km domain shows little change in 
stagnation in the northeast, the high resolution simulation shows stagnation increases 
of up to 5 days per season during the summer over parts of the northeast which is large 
compared to the average 5 to 10 stagnation days per season that occur in this region 
20 (Fig. 5b). During fall, a large increase (over 1 5 days) in the number of stagnation days is 
found along the Gulf coast and the California coast. Along the coast of California during 
fall, the increase in stagnation days leads to increased concentrations of pollutants, 
reinforcing the negative impact that increased temperature and insolation have on air 
quality in the area. Similarly, the decreased precipitation along the Gulf coast may 
25 reinforce higher concentrations of pollutants due to increased stagnation in the area. 
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The ventilation coefficient, which is defined as the product of the mean wind speed 
within the boundary layer and the boundary layer height (Pielke et al., 1991), reflects 
how well pollutants can be mixed and transported within the boundary layer. Ventila- 
tion is adversely impacted by stagnation, which is driven by the persistence of certain 
large-scale circulation patters, but also takes into account smaller scale meteorological 
conditions. An unvented hour is an hour during which the ventilation coefficient is less 
than 6000m 2 s _1 . During summer, Texas has on average 1 to 2 more unvented hours 
per day in the future compaired to the present which tends to increase pollutant concen- 
trations, further amplifying the increased concentrations of 0 3 and some secondary PM 
(with the exception of volatile PM such as ammonium nitrate) due to increased temper- 
atures in the region (Fig. 6a). The 12 km simulation shows that the coast of Georgia and 
the Carolinas also see 1 to 2 more unvented hours per day during the summer, while 
the 36 km shows less dramatic changes over the southeast. During fall, an increase of 
1 to 2 unvented hours per day is found over much of the northeast and southeast. The 
higher resolution domains show similar trends, although with more spatial variability. In 
the northeast, the combined higher temperatures and less ventilation would influence 
higher concentration of pollutants, while the decrease in insolation would reduce sec- 
ondary pollutants such as 0 3 and secondary PM (Fig. 6b). Unvented hours over most 
of Minnesota increase during spring by over 3 h per day; however, none of the other 
variables examined here show either a reinforcing or competing effect on air quality. 

4.5 Regional climate and urban centers 

U)' 

Since a large and growing fraction of the nation’s population is located in dense urban 
areas, it is important to examine the change in air quality related climate variables over 
some major US cities. The expected response to climate change differs among var- 
ious different regions of the US. Here we focus on 5 geographically unique, densely 
populated cities that are representative of the different regions of the US: Atlanta, 
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Los Angeles, Philadelphia, Phoenix, and Seattle. Further, land use changes may ex- 
acerbate meteorological changes in cities. Philadelphia and Phoenix are chosen here, 
rather than more populated cities such as New York and Houston, because future stud- 
ies are planned to address the impact of land use changes on regional climate in those 
5 cities. Extremes in meteorological variables are important because there are adverse 
health effects associated with short term exposure to poor air quality. Air pollution is 
highly variable in time and temperature extremes are also associated with adverse 
health outcomes (McMichael et al., 2006). Cumulative distribution function (CDF) plots 
show the percentage of hourly temperature and daily precipitation that exceed a given 
10 value in major US cities for each simulation year (Fig. 7a-e). The hourly temperatures 
within any given percentile range tend to shift a few degrees warmer in the future, ex- 
cept for the lower 20th percentile range in Atlanta and Philadelphia and the upper 90th 
percentile range in Seattle. Most warming in Seattle occurs at the lower 75 percentile 
range, where high 0 3 concentrations are not likely, which reflects the increase in tem- 
i5 perature mentioned earlier during the winter in the Pacific Northwest. This can also 
decrease emissions related to domestic heating, including PM from wood burning. The 
cumulative distribution of the maximum daily temperature in Seattle follows a similar 
trend, as the hourly CDF in the lower range and temperature are nearly the same in 
the 60 to 90 percentile range (Fig. 8a-e). The upper 95th percentile in daily maximum 
20 temperatures in Seattle are actually around a degree cooler in the future, decreasing 
the chance of high ozone during hot days. Los Angeles and Phoenix have similar hourly 
temperature and maximum temperature CDF structures and shifts from present to fu- 
ture. These cities are warmed by 1 to 3 degrees regardless of the percentile range. Lin 
et al. (2001) have developed estimates of the probability that the maximum daily 8h 
25 average 0 3 will exceed 80ppb given the maximum daily temperature in a given region 
(including Los Angeles, the southeast, and the northeast). Given that the upper 95th 
percentile in daily maximum temperatures in Phoenix increases from around 308 K to 
31 2 K, the probability that 0 3 will exceed 80ppb on these days increases. Similarly, 
a shift in the upper 95th percentile in daily maximum temperatures in Los Angeles, 
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from around 305 K to 308 K, also increases the probability of high 0 3 . The coolest days 
in Atlanta and Philadelphia are similar to, if not cooler than the future simulations, re- 
flecting the cooling that occurs during the winter and spring in the surrounding regions. 
While very little, if any, warming occurs in the lower 50th percentile of daily maximum 
5 temperature, the upper 95th percentile maximum temperatures in Philadelphia and 
Atlanta increase from around 302 K to 303 K and 304 K to about 307 K respectively, 
with a corresponding increase in probability of high 0 3 concentrations on those hot 
days. The 12 km simulations produce similar, but not identical, cumulative distribution 
of total hourly temperatures and maximum temperatures in Atlanta and Philadelphia 
io (Fig. 9). The 12 km domains show a shift in daily maximum temperature at the upper 
95th percentile of 303 to 305 K and 304 to 305 K in Philadelphia and Atlanta respec- 
tively, implying a higher probability of high ozone on those days than would be derived 
from the coarse resolution. The hottest days in Philadelphia are simulated by the fu- 
ture 12 km domain where temperatures reach 310 K. Maximum temperatures and high 
is 0 3 probability in New York exhibit similar changes to that of Philadelphia. Temperature 
distributions in Chicago were also analyzed, however the only difference from present 
to future is a shift in the upper 95th percentile from 300 to 301 K while the distribution 
below the 90th percentile does not change. 

Rain can improve air quality, and although the seasonal mean precipitation has al- 
20 ready been examined, it is important to also understand how the frequency and distri- 
bution of rainfall can shift over time. Rainfall frequency shifts are especially important in 
cities, where frequent light rains will clean the air more than infrequent heavy rains. The 
most notable change in daily precipitation distribution occurs in Los Angeles, where the 
driest year is simulated in the future and the wettest year is simulated in the historic 
25 simulations (Fig. 10). Seattle does not seem to receive more or less rain in the fu- 
ture but does have a more consistent distribution of daily rainfall from year to year in 
the future. The precipitation distributions of the remaining cities do not appear to be af- 
fected. Precipitation distributions produced from the 12 km simulations over Atlanta and 
Philadelphia are similar to the distribution produced by the coarse simulation (Fig. S2). 
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5 Conclusions 

O 

C 

The high resolution simulations produce somewhat different results than the coarse 
resolution simulations in some regions. Also, through the analysis of the meteorological 
variables that most strongly influence air quality, we find consistent changes in regional 
climate that would enhance ozone levels in four regions of the US during fall, one 
region during summer, and one region where changes potentially would lead to better 
air quality during spring. 

During summer and fall, all air quality indicators, with the exception of insolation and 
precipitation, suggest an increase in air pollutant concentrations, including increased 
production of secondary PM and 0 3 , in most of Texas. Consistent with Leung and 
Gustafson (Leung and Gustafson, 2005), when comparing 2048-2052 to 2006-2010, 
we find warmer temperatures (2-3 K), less ventilation (1-2 h per day) and more stag- 
nation (10-15 days per season) during summer in Texas and the same for fall, though 
slightly less intense. 

The West Coast is marked by warmer temperatures (ranging from 1-4 K), slightly 
less rainfall (less than 1 mm per day), and more stagnation (10-15 days) during fall, 
while there is no significant change in ventilation and insolation. Similarly during fall, the 
southeast shows little change in ventilation and stagnation but is warmer (2-3 K), dryer 
(up to 2 mm per day in some areas) and with slightly higher insolation (5-1 OWrrf 2 ). 

Also, the shift in the distribution of maximum daily temperatures in Atlanta increases 
the probability of high 0 3 concentrations for days when the maximum temperature is in 
the 95th percentile. 

The northeast is also warmer during fall (2-3 K) and less ventilated (1-2 h per day) 
but increased rainfall (up to 2 mm per day) and decreased insolation (~5Wrrf 2 ) com- 
pete for better air quality. Regardless, as in Atlanta, the shift in the distribution of maxi- 
mum daily temperature in New York and Philadelphia increases the probability of high 
0 3 concentrations (0.06-0.12) for days when the maximum temperature is in the 95th 
percentile. During spring, however, increased rain (1-2 mm per day), more ventilation 
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(~ 1 h per day) and decreased temperatures (~ 1 K) could promote better air quality in 
the northeast. 

< / > 

While climate conditions strongly impact air quality, emissions and chemistry also 
play a vital and complex role in the formation and removal of atmospheric pollutants. 

5 A more comprehensive assessment of emissions and chemistry will be addressed in 
the future. 


Supplementary material related to this article is available online at: 
http://www.geosci-model-dev-discuss. net/6/251 7/201 3/ 
gmdd-6-251 7-201 3-supplement.pdf. 
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Table 1. GISS-WRF modeling system performance for wind speed and direction against TDL 
observations for US regions and seasons. 
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Table 2. GISS-WRF modeling system performance for temperature against TDL observations 
for US regions and seasons. 
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Fig. 1 . Modelling domains with horiontal grid-spacing resolutions of 36 km and 1 2 km northeast 
(NE) and southeast (SE). 
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Fig. 2. Predicted average yearly and seasonal 2 m atmospheric temperature change (future 
minus historic) for (a) the 36 x 36km resolution modeling domain, (b) the 12 x 12km resolution 
sub-domain over northeast and (c) the 12 x 12 km resolution sub-domain over southeast. 
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Fig. 3. (a) Predicted average yearly and seasonal precipitation (mm day -1 ) change (future 
minus historic) for (a) the 36 x 36km resolution modeling domain, (b) the 12 x 12km resolution 
sub-domain over northeast and (c) the 12 x 12 km resolution sub-domain over southeast. 
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Fig. 4. Predicted average yearly and seasonal downward short wave radiative flux at the surface 
(W rrf 2 ) change (future minus historic) for(a)the 36x36 km resolution modeling domain, (b)the 
12 x 12km resolution sub-domain over northeast and (c) the 12 x 12km resolution sub-domain 
over southeast. 
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Days 


Fig. 5. Predicted total seasonal change in the number of stagnation days (days per season) 
(future minus historic) for (a) the 36 x 36km resolution modeling domain, (b) the 12 x 12 km 
resolution sub-domain over northeast and (c) the 12 x 12km resolution sub-domain over south- 
east. 
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Fig. 6. Predicted total seasonal change in average unvented hours (hours per day) (future 
minus historic) for (a) the 36 x 36km resolution modeling domain, (b) the 12 x 12km resolution 
sub-domain over northeast and (c) the 12 x 12 km resolution sub-domain over southeast. 
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Fig. 7. Empirically determined cumulative distribution of 36 km historic (2006-2010) and future 
(2048-2050) hourly temperatures at major US cities: (a) Atlanta, (b) Los Angeles, (c) Philadel- 
phia, (d) Phoenix and (e) Seattle. 
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Fig. 8. Empirically determined cumulative distribution of 36 km historic (2006-2010) and future 
(2048-2050) maximum daily 1 h average temperature at major US cities: (a) Atlanta, (b) Los 
Angeles, (c) Philadelphia, (d) Phoenix and (e) Seattle. 
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Fig. 9. Empirically determined cumulative distribution of 12 km historic (2006-2010) and future 
(2048-2050) hourly temperature (a) Atlanta, (b) Philadelphia and maximum daily 1 h average 
temperature at (c) Atlanta, (d) Philadelphia. 
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Fig. 10. Cumulative distribution of 36 km historic (dark) and future (light) daily precipitation at 
major US cities: (a) Atlanta, (b) Los Angeles, (c) Philadelphia, (d) Phoenix and (e) Seattle. 
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